Theory of Scanning Tunneling Spectroscopy of Magnetic Adatoms in Graphene 
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We examine theoretically the signatures of magnetic adatoms in graphene probed by scanning 
tunneling spectroscopy (STS). When the adatom hybridizes equally with the two graphene sublat- 
tices, the broadening of the local adatom level is anomalous and can scale with the cube of the 
energy. In contrast to ordinary metal surfaces, the adatom local moment can be suppressed by the 
proximity of the probing scanning tip. We propose that the dependence of the tunneling conduc- 
tance on the distance between the tip and the adatom can provide a clear signature for the presence 
of local magnetic moments. We also show that tunneling conductance can distinguish whether the 
adatom is located on top of a carbon atom or in the center of a honeycomb hexagon. 

PACS numbers: 73.20.Hb,71.55.Ht,73.20.-r 
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Graphene is a two dimensional sheet of carbons whose 
remarkable electronic properties derive from the pres- 
ence of electronic excitations that behave as chiral Dirac 
quasiparticles[l]. Although clean bulk graphene may not 
be magnetic, there is a rich variety of possibilities for 
magnetism when adatoms are added on top of graphene. 
As an open surface, the use of scanning tunneling mi- 
croscopy (STM) probes [2[ opens the possibility of con- 
trolling the position of adatoms with atomic precision [3] 
and at the same time switching the magnetic local mo- 
ments on and off by gating 0, 

One of the challenges for the manipulation of local mo- 
ments is about detection: how can one reliably identify 
a local magnetic moment at room temperature. Unlike 
ordinary metal surfaces, due to the low density of states 
(DOS), an adatom localized level can hybridize strongly 
with the STM tip. We propose that the dependence 
of the adatom STM differential conductance (DC) with 
the distance to a non-magnetic STM tip, in the limit 
of small separation, can provide an experimental signa- 
ture for the presence of local moments above the Kondo 
temperature jfj. Furthermore, because the electrons in 
graphene carry different sublattice quantum numbers, we 
show that when the adatom sits in the center of the hon- 
eycomb hexagon [see Fig. 1(b)], destructive interference 
between the different tunneling paths changes substan- 
tially the form of the Fano factor [7[ and the shape of the 
DC curves compared to the case where the adatom is lo- 
cated on the top of a carbon atom [Fig. 1(a)]. This effect 
allows the use of STM to characterize adatoms and de- 
fects in graphene, including substitutional impurities in 
single and double vacancies [8]. 

Our starting point is the free Hamiltonian of the 
graphene/adatom/tip system: TL = H g +Hf+H c - H g 
is the tight-binding Hamiltonian for graphene: TL g = 
Sa a i( R & o-(Rj) + h.c, where a, b are the 
fermionic operators for sublattices A and £>, respectively 



(t ~ 2.8 eV), and a =jj is the spin. In momentum space, 



n 9 = -tj2 kk)4A. + f(k)i« kfl 



(i) 



where 0(k) = E;=i a i = 8l 2 = -x/2 + V^y/2, 

and 3-s = —x/2 — \/3y/2 are the lattice nearest neighbor 
vectors. TL C = J2k e k<4- ^Cka is the effective Hamilto- 
nian for the c tip electrons, with ek = (/c 2 )/2m* — ep, 
where m* is the effective electronic mass, and ejj is the 
energy at the bottom of the tip band with respect to the 
Dirac point, and Hf = to f a fo- is the Hamiltonian 
of the / electrons at the local level with energy eo- The 
Coulomb energy, U, for double occupancy of the local 
level is described by a Hubbard term: TLu = C^/|/t/|/|- 
Since we are only interested in the magnetic state of the 
adatom (we do not include the Kondo effect and hence 
our theory is valid above the adatom Kondo tempera- 
ture, Tr-), in what follows we use Anderson's mean field 
decomposition [4]: Hu,mf = U Ysa n °f-af-o- ~ Un^, 
where n a = (f£f a ) is the average occupation of the 
level. Hence, at the mean field level we write: 7Y/,mf = 
X^cr e a fa fa ? where e a — eo + Un- a is the renormalized 
level energy. 

In graphene the adatoms can be localized at different 





FIG. 1: Two adatom positions in graphene: a) asymmetric 
case, on top of a carbon atom, when the adatom (large circle) 
hybridizes with one sublattice and b) symmetric case, when 
the adatom is located at the center of a hexagon. Red arrows: 
nearest neighbor vectors. 
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positions in the honeycomb lattice [8[. Here we consider 
two cases where the adatom is either placed on top of 
a carbon atom, in which case the sublattice symmetry 
is locally broken, or the adatom is located in the center 
of an hexagon without symmetry breaking. In the first 
case, assuming the adatom to be on the B sublattice, we 
have: 7Yy,AS = V f^b a (0) + h.c. In the second case 

we have: H v ,s = J2aYfi=i V [ a a( a + & a(~ a 0] f* + h.c. 
In momentum space, these two terms can be written as: 

Wv,AS=V^6t (7 / <T + h.c. (2) 

per 

in the asymmetric case, and 

H v ,s = V^ W)*U + U + h-c (3) 

per 

in the symmetric one [q|. 

In the presence of an STM tip, there are two additional 
hopping terms: Hf- C = V c ]T a fa c <j( r ) + n - c -> where V c 
is the effective tunneling energy between the tip and the 
adatom, and H g - C = Y, a ,i *c,i [^(Ri) + Ca(Ri- 
r) + h.c. where t C) i is the hopping energy between the tip 
and each carbon, r = (i?, z) is the position of the tip with 
respect to the adatom, R is the in-plane distance from the 
adatom to the center of the tip, and z is the out-of-plane 
distance between them. The fact that the adatom sits a 
few angstroms above the graphene plane is accounted by 
the exponential z dependence of the tunneling energies. 
In momentum space we have: Hf- C = V c (r) ^2 ap f^Co-p J r 
h.c, and H g - C = Eakp^pO) [4 P + ^ P ] c ak + h.c, 
where t CjP (r) = t c (z)e~ ip '' R after averaging over the po- 
sition of the carbon sites below the impurity. 

In the absence of the tip, after diagonalization of the 
Hamiltonian H g + Hf + H v in Eq. Q-©, the re- 
tarded Green's function of the / electrons Gff a (r) = 
-(T[/(r)/t(0))])isgivenby: 

Gf f>< ,{u) = [u-e*-'Eff(u)r 1 (4) 

where Y>ff(uj) is the self-energy of the localized electrons 
due to the hybridization with the electrons in graphene. 
As in the usual Anderson impurity problem[10], the 
formation of local moments is defined by the occupa- 
tion of the localized level n a for up and down spins, 
n a — — (1/tt) f^ D duImGf-f a (uj), where \i is the chemi- 
cal potential (D ~ 7 eV is an energy cut-off). 

In the asymmetrical case [see Eq.flS])], Ejj(cj) = 
V 2 E p G° bb R p (u;), where G° bb (r) = - (T[b(r)tf (0)}) is the 
diagonal component of the bare graphene Green's func- 
tion, G^y, with x, y = a, b: 

r 0R , x = U(J ly ~ ^yRe0(k) + ftr^Im^k) 
^' klCJj u 2 - t 2 |(/>(k)| 2 + zO+sign(a;) ' 1 } 

where a J (j = 1,2) are off-diagonal Pauli matrices (a 2 a = 
a ab ~ 0? anc ^ a ° 1S the identity matrix. We calculate the 



self-energy within the linearized theory, where the energy 
around the K (K' ) point is given by: e±(k) = ±vpk, 
where vp « 10 6 m/s is the Fermi velocity. It reads [4]: 

E^M = —(vA/tt) In (\D 2 - co 2 \/co 2 ) -iA\co\0(D , 

(6) 

where A = tt(V/D) 2 is the dimensionless hybridization 
parameter. The imaginary part of (|6|) gives the broaden- 
ing of the level due to the hybridization of the adatom, 
and it is proportional to the DOS of the host. 

In the symmetric case, however, the hybridization of 
the / electrons is mediated by the virtual hopping of the 
electrons in graphene into a "ghost" site located below 
the adatom, in the center of the hexagon. In this process, 
as the electrons hop in and out from the adatom, they 
gain an additional phase that leads to interference be- 
tween the different quantum mechanical paths involving 
the six adatom neighboring sites on both sublattices. The 
self-energy in this case involves also off-diagonal terms of 

e^M = v 2 J2 Wp)e 0)P M + 0*(p)e b , P H] , (7) 

P 

where e x , p (u>) = cj>(p)G xb R p (cj) + f(p)G x R p (u). In the 
linearized theory, Eq. ([7j) gives 11 : 

= -u [Z' 1 ^) - 1] - 2iA(\co\ 3 /t 2 )0(D - H) , 

(8) 

where Z' 1 = 1 - (l/wjReE^w) = 1 + 2A/(irt 2 )[D 2 + 
u 2 In [\D 2 — u 2 \/u 2 ^} gives the quasiparticle residue. 
The imaginary part of Ejj gives rise to an anomalous 
broadening of the adatom level that scales with |cj| 3 /t 2 , 
suppressing strongly the hybridization when |eo| <t. 

In the perturbative regime where V c , t c are small com- 
pared to V, the inclusion of the tip leads to an addi- 
tional renormalization of the / electrons Green's func- 
tion, £//(r, uS) = Tiff{uS) + Eyrj (r, a;), where 

E«(r, w ) = V c (r,u)V c (-r,co)J2G° cc %(u;). (9) 

k 

G ccfkM = - e(k) +i0 + ] _1 is the retarded Green's 
function of the c-electrons, G° cc1l (t) = — (T[ck(r)cj c (0)]), 
whereas V c (r, uS) is the renormalized tunneling energy be- 
tween the tip and the adatom, namely 

V c as (t,oj) = V c + Vt c (z)[G a R (K,w)+G° b R (R,u;)] (10) 

for the asymmetric case, where V c (r, uj) follows from the 
exchange G® xh — > G° x , and 

V c s (r, u)=V c + Vt c (z) [6 a (R, lj) + 9 6 (R, u)] (11) 

in the symmetric one (V s = V s ). In our notation, 
6(R) = Xlp e * p R ®^,P i s the Fourier transform of 6^. 
Ey-j in Eq. © can be easily computed assuming an ef- 
fective band width, ap, for the c electrons in the tip. 



3 




0.2 0.4 0.6 

bias (eV) 



0.214 0.216 0.218 

bias (eV) 



FIG. 2: Adatom induced DC, Q =dI/dV, versus bias when 
the adatom sits (left) on top of carbon and (right) in the 
center of the hexagon. See details in the text, a), b) t c = 
0.15eV and V c /t c = 1.6, 1.5, 1.35, 1.1, 0.7, 0.25 and 0.1 (inset), 
from top to bottom, c), d) t c = 0.02eV and V c /t c = 0.7, 0.25, 
0.1, 0.05 and 0.01. V c /t c = 0.1 for all curves in red. 



The adatom induced DC is £(r, uo) oc p c (uo)pf- c (r, uo), 
where p c is the DOS of the tip, and p/_ c (r, uo) = 
-(l/7r)Im \V c (y,uj) ^ a G f f^(uo)V c (-Y,uo)] is the ana- 
logue of the / electron DOS, which contains renormal- 
ized tunneling matrix elements between the tip and the 
adatom [l2[. In a more standard form, the DC is 



Q(t,cj) = 2irep i 



,2 ~ QQ 



■ 77 + (<27 + Ql)£o 



£ 2 



l 



(12) 



= \i + eV, with 
2 (u)/lmV ff (r,w), 



where e is the electron charge, uo 
eV the applied bias, p(r,u) = —p 
where p(uo) is the graphene DOS per spin, and 
€a(r,uj) = -[uo - e a - Re£//(r,cj)]/Im£//(r,cj), 
q(r,uo) — ReV r c (r, uo)/[7rt c (z)Vp(uo)] and 7(1*, uo) = 
— ImT4(r, uo)/[7rt c (z)Vp(uo)] are the Fano parameters Q. 
In contrast with usual metal surfaces, the Fano factor q 
has a live dependence with the bias in graphene. The 
conjugate forms q(— r, uo) and j(— r, uo) are defined by 
the conjugate tunneling matrix element V c (— r, uo). In 
the particular case when the tip is above the adatom 
(R = 0), q = q and 7 = 7 and q AS = [V c + 
(t c (z)/V)ReEff(uo)}/[7Tt c (z)Vp(uo)}, where is 
defined by Eq.(j6]), and j AS = 1. In the case where the 
adatom and the tip are on top of each other in the center 
of the hexagon, destructive interference leads to cancel- 
lation of the perturbative corrections on the tunneling 
matrix element in Eq.([TT]h and the Fano parameter sim- 
plifies to q s = V c / [7rt c (z)Vp] and 7^ = 0. For adatoms 
with d and /-wave orbit als, the cancellation is not exact. 

The shape of the Fano resonances in the DC curves 
is driven by the ratio q/j. When q/j ^> 1, the DC 



curve shows a pronounced peak, whereas in the opposite 
regime, q/j <C 1 one expects a dip. For a set of pa- 
rameters V = leV, U = leV, ajj = 4 eV, ejj = 2eV, 
p = O.leV and eo = — 0.5eV, in the asymmetric case, for 
Vc/tc = 0.1, the red curve shown in the inset of Fig. 2(a) 
has a small dip, which is suppressed when V c /t c > 0.2. 
In contrast, the curves shown in Fig. 2(b) for the sym- 
metric case have a well pronounced peak for all finite 
values of q s [see inset of Fig. 2(b)], reflecting the fact 
that q s /j s is always large (j s = 0). Fig. 2(a),(b) and 
2(c),(d) compare the features of the positive bias reso- 
nance for t c = 0.15eV and t c = 0.02 eV, respectively. All 
red curves in Fig. 2 correspond to V c /t c = 0.1. Increasing 
this ratio from V c /t c = 0.25 up to 1.6, the DC curves 
show two strongly pronounced peaks indicating the po- 
sition of the two magnetic Fano resonances at eo + n^U 
and eo + n^U. For V c /t c < 0.1, in the asymmetric case, 
Fig. 2(c) shows an inversion in the structure of the res- 
onance (Vc/tc = 0.01) for positive bias. The shape of 
the DC curves for small q/j [see insets of Fig. 2] agrees 
with a recent STM measurement of the Kondo peak in 
graphene [l^. 

The decrease in the separation of the peaks with in- 
creasing V c reflects the suppression of the local moment 
by the proximity of the STM tip. In particular, in the 
symmetric case [Fig. 2(b)], the hybridization of the 
adatom with graphene is weaker than in the asymmetric 
one, making local moment much more sensitive to the 
STM tip. The difference is indicated clearly by the sepa- 
ration of the peaks for large V c and also by the width of 
the peaks as V c goes to zero. In the asymmetric case (Fig. 
2(a)), the peaks remain broad at small V c whereas in the 
other case their width collapses much faster, reflecting 
their anomalous broadening oc \uo \ 3 /t 2 [see Eq. ©]. In 
the opposite limit, for V c large enough, the two DC mag- 
netic peaks eventually merge on top of each other, de- 
stroying the local moment completely. The merging of 
the peaks happens much earlier in the symmetric case 
[Fig. 2(b)] than in the asymmetric one. We note that for 
large V c , fluctuations drive the DC away from equilibrium 
(although still in the perturbative regime for V c <C V), 
invalidating the strict applicability of Eq. ([T2]) . The 
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FIG. 3: Graphene LDOS at the adatom site (top carbon case) 
for Vc/t c = 0.1, 0.7, 1.1 and 1.6, from (a) to (d) (t c = 0.15eV). 
Black curve: n^; brown: nj.. Total LDOS: n-f +7i|. 
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FIG. 4: Integrated LDOS maps around the adatom (center) 
in the asymmetric case, with the adatom on top of the carbon. 
On the left (right): scans for the opposite (same) sublattice 
of the adatom: (a),(c) without (t c = V c = 0) and (b),(d) with 
the tip (t c = 0.2eV, V c = 0.02eV). 



main effect of fluctuations, however, is to further sup- 
press the local moment, and in this sense the equilib- 
rium calculation may be regarded as a conservative es- 
timate for the main effect, the suppression of the local 
moment by the STM tip. As the tip separation to the 
adatom becomes progressively small the DC peaks can 
shift strongly: in that case, the peak on the right will 
red shift and eventually cross the experimentally acces- 
sible bias window around the Fermi level, providing an 
experimental signature for the presence of the local mo- 
ment with a non-magnetic tip, regardless the presence of 
the Kondo peak. A significant suppression of the local 
moment by the metallic tip is harder for adatoms with 
a very large U, such as cobalt, which show a large local 
moment when hybridized with metals, but may be easily 
achieved in adatoms which are not usually magnetic and 
exhibit a local moment in graphene [4|. 

The Fano resonance also generates magnetic peaks in 
the LDOS around the adatom. The plots in Fig. 3 show 
the evolution of those peaks with the increase of V c in 
the asymmetric case, for up and down spins. As the level 
is broadened by the tip, the height of the peaks clearly 
collapses [Fig. 3(d)]. 

In the case where the adatom is on top of the carbon, 
the sublattice asymmetry in the LDOS provides another 
STM signature that identifies the position of the adatom 
in the lattice. The LDOS for the asymmetric case can be 
computed perturbatively for small t c and V c , 

(£S(r,u>) = -l/(KV 2 )ImZf^(z,Lj)\ Vc=0 -(lMx 
Im [T x (r,w)Gf fia (w)f x (-r,Lj)] , (13) 

where x = b for the same sublattice of the adatom and 
x = a for the opposite one. The first term is the renor- 



malized LDOS in the absence of the adatom and 

I x (r, u) - VG xb (K, u) + v^ s {-t,lj) ff ^ ' 

contains the interference effects due to the interplay of 
the adatom and the tip in graphene. T x in Eq. ([T3]) fol- 
lows by exchanging G xb by G bx and V AS by V AS [see 
Eq. ([T0|) ]. The STM topography maps computed from in- 
tegration of Eq. (p~3|) in energy (see Fig. 4) clearly show 
the asymmetry between the two sublattices (the adatom 
has three nearest neighbors and six next nearest neigh- 
bors). Fig. 4(a), (b) display the integrated LDOS for 
the opposite sublattice of the adatom, which has a lower 
point group symmetry, while Fig. 4(c), (d) display similar 
maps for the same sublattice of the adatom, with and 
without interference effects from the tip. 

In conclusion, we have derived the fingerprints for Fano 
resonances of magnetic adatoms in graphene. We have 
shown the signatures in the DC curves that identify the 
position of the adatom and possibly the presence of local 
moments, away from the Kondo regime. 
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